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Abstract 

A theory for longitudinal (Ti) and transverse (T2) electron spin coherence 
times in zincblende semiconductor quantum wells is developed based on a 
non-perturbative nanostructure model solved in a fourteen-band restricted 
basis set. Distinctly different dependences of coherence times on mobility, 
quantization energy, and temperature are found from previous calculations. 
Quantitative agreement between our calculations and measurements is found 
for GaAs/AlGaAs, hiGaAs/InP, and GaSb/AlSb quantum wells. 
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The recent observation of long-lived (> 100 ns) spatially extended (> 100 fim) coher- 
ent spin states in semiconductors suggests the possibility of manipulating nonequilibrium 
electron coherence to an unprecedented degree in a solid, [jl]-^] These spin states interact 
with light, and thus can be used to generate a host of novel dynamic nonlinear optical and 
electrical effects. The magnitude and persistence of such effects is governed partly by 
the spin coherence times Ti and T2, describing the decay of longitudinal and transverse spin 
order, respectively. Ultrafast optical measurements have been performed of both Ti and T2, 
although in different geometries p|-^,|6|-[T0[ . 



To guide further efforts in the controllable manipulation of room-temperature spin co- 
herence it is essential to have a quantitative theory of spin decoherence. Direct quantita- 
tive comparison of the current theory with experiment has been rare for quantum wells, 
for an independent measurement of the mobility of the quantum well is required. Re- 
cently such a comparison was made for room-temperature electron spin lifetimes in n-doped 
GaAs/AlGaAs multiple quantum wells (MQWs) 0. In addition to measured Ti's one order 
of magnitude longer than those predicted from current theory there were discrepancies in 
the power law dependences of Ti on mobility and confinement energy. 

This Letter provides the desired quantitatively accurate theory of spin decoherence for 
quantum wells and clarifies the relationship between Ti and T2 in these systems. Our results 
are in excellent agreement with experimental measurements on GaAs quantum wells not 
only in the previously unexplained general trends, but also in the absolute magnitude. We 
also find excellent agreement with measurements on InGaAs/InP and GaSb/AlSb 
quantum wells, whereas previous calculations disagree by an order of magnitude. Finally 
we find unexpected trends in the spin coherence times with temperature which may explain 
other puzzling experimental results. 

The mechanism of electron spin decoherence we consider occurs via the spin precession of 
carriers with finite crystal momentum k in the effective k-dependent crystal magnetic field 
of an inversion-asymmetric material. A signature of this mechanism is that in the "motional 
narrowing" regime, where orbital scattering times r greatly exceed spin decoherence times. 



Ti oc r~ . Thus cleaner samples have shorter spin coherence times. In III-V bulk ||TT 
and quantum well structures nanostructures this trend has been observed in samples of 
varying mobility at and near room temperature. 

D'yakonov and Perel' have developed a theory for Ti based on this mechanism for bulk 
zincblende semiconductors, assuming orbital coherence is lost after each scattering event, 
and assuming a thermal distribution of electrons |T2|. This work was later extended (with 
further approximations) to quantum wells by D'yakonov and Kachorovskii (DK). Thus 
the two categories of approximation in DK theory are (1) the method of handling the 
orbital degrees of freedom and (2) the quantum well electronic structure. For example, 
if some orbital coherence or nonthermal occupation were maintained after each scattering 
event, then because the electron's orbital degrees of freedom are entangled with its spin, 
only a nonequilibrium calculation of orbital degrees of freedom (e.g. Monte Carlo) would 
produce accurate results. We find, however, that the sources of error are the approximations 
pertaining to quantum well electronic structure, and thus simpler, quantitatively accurate 
calculations of spin coherence may be performed. 

Our theory begins with the assumption of motional narrowing. In the motional narrowing 
regime the electronic spin system is subject to an effective time- dependent, randomly oriented 
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magnetic field H which changes direction with a time r that is much shorter than the 
precession time of either the constant apphed field Ho or the random field. The coherence 
times depend on the transverse {Hj_) and longitudinal {H\\) components of the random field, 
according to fl^ 

Tfi oc {Hl)T, (1) 
T,-V([i7i]/2 + i7f)r, (2) 

where the constant of proportionality is the same for Eqs. (|l]) and (^. In a crystal with 
inversion asymmetry and spin-orbit coupling there is a spin splitting described by the Hamil- 
tonian H = hft(k) ■ a/2, where f2(k) is a momentum- dependent eSective magnetic field. As 
the electron is scattered from k to k' via ordinary orbital (not spin- dependent) elastic scat- 
tering, the effective magnetic field changes direction with time. If the crystal is cubic, then 
= Hy = iff, so Hj^ = 2H^^ and T2 = Ti. The relationship between Ti and T2 differs, 
however, for systems of lower symmetry, such as quantum wells. For a (001) grown quantum 
well the fluctuating field along the growth direction vanishes, and 

Tf 1(a) = T^\a = 0)(1 + cos^ a)/2, (3) 
T^\a) = T^\a = 0)(2 + sm^a)/4, (4) 

where a is the direction between Ho and the growth direction. Thus T2 ranges from 2Ti/3 
to 2Ti depending on a. In contrast, for (110) grown quantum wells the effective crystal 
magnetic field is entirely along the growth direction, and 

Tf 1(a) = T^\a = 0) sin^ a, (5) 
T^\a) = T^\a = 0)(1 + cos^ a)/2. (6) 

Thus although T{~^{a = 0) vanishes, the same is not the case for Tg"! for any a. 

Calculations of Ti require knowledge of both r and f2(k). The effective time for field 
reversal (r^) depends on the angular index i of the field component (fi^). For example, an 
i = 1 component {fli) requires a 180° change in the angle of k to change sign, whereas an 
i = 3 component (^3) only requires a 60° change, so typically < ti. Thus 

^ = - f D{E)f{E)[l~f{E)]Y^n{E)n]{E)dE, (7) 

where f{E) is the Fermi occupation function, D{E) is the density of states, n is the elec- 
tron density, and the scattering rates t^^{E) = f\a{9, E){1 — P£{cos6))dcos9 for bulk 
semiconductors and t^^{E) = Q'^ a{9,E){l — cos[l9])d9 for (001) quantum wells. For both 
bulk and quantum wells the functional form of the scattering cross-section a{9, E) is taken 
from standard expressions for ionized impurity (II), neutral impurity (NI — such as arises 
from quantum well interface roughness), or optical phonon (OP) scattering. The r^'s differ 
for different mechanisms (e.g., for a quantum well Ti/t^ = i"^ for II, Ti/re = i for OP, and 
Ti/Tf = 1 for NI scattering). The magnitude of a{9,E) is obtained from the mobility, 

/i = (e/mn) J D{E)f{E)[l - f{E)]n{E)EdE. (8) 
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We obtain fle{E) from a non-perturbative calculation in a fourteen-band basis [|T5[. This 



basis, which is the minimum required to generate spin splitting nonperturbatively, consists 
of two conduction antibonding s states (s), six valence (bonding) p states, and six anti- 
bonding p states (p). Such a basis has, for example, been used to analyze spin-splitting in 
heterostructures ||T6|. The Hamiltonian is well-known, and can be found in Refs. |T^,|T5[. 
The parameters that enter this Hamiltonian include the zone-center energies of the con- 
stituent bulk semiconductors and the momentum matrix elements between bands, which 
are obtained from the conduction band mass, the heavy-hole mass, and the (7-factor. Time 
reversal invariance requires ri(k) = — ri(— k), so = for even £. 

For quantum wells the electronic structure is obtained by expressing the electronic states 
as spatially-dependent linear combinations of the fourteen states in the basis. The full Hamil- 
tonian is projected onto this restricted basis set, which produces a set of fourteen coupled 
differential equations for the spatially-dependent coefficients of the basis states (generalized 
envelope functions). These equations are then solved in Fourier space in a similar method 



to that of Winkler and Rossler. [|18| Further details are available elsewhere [jlS . 

For bulk semiconductors the relevant electronic states for spin decoherence are near the 
bulk band edge, and thus perturbative expressions for fl'j{E) for these bulk semiconductors 
(QKE) = 0, f^K-E) oc E^) |[T^ are identical to those obtained from a full fourteen-band 
calculation within numerical accuracy. Shown in Fig. |l| are calculated T2's for GaAs, InAs, 
and GaSb assuming II scattering. The agreement with experimental measurements for 
GaAs at the higher temperatures is quite good, whereas for low temperatures other spin 
relaxation mechanisms are expected to dominate. The smaller T2's in InAs and GaSb are 
due partly to the larger conduction spin splitting, which originates from a larger ratio of the 



spin-orbit coupling A to the band gap Eg (see Ref. |]T9| on perturbative expansions of spin 
splittings). The agreement between calculated and measured T2's in Fig. |I] indicates that 
the spin splitting of bulk GaAs is well described by our model. 

We now contrast our results for quantum wells with those of DK theory. The DK theory 
for (001) quantum wells is derived as follows. First, negligible penetration of the electronic 
states into the barriers is assumed, so Ei -C AEc, where Ei is the confinement energy of 
the first quantum well state and AE^ is the conduction band offset. Then the perturbative 
expressions pO|] ^1{E) oc E{4Ei—E)'^ and ^"^{E) oc E^ are used. Furthermore the energies of 
relevant states are assumed to be <C -Ei, and thus (i) the contribution from Qsi^E) is ignored, 
and (ii) it is assumed that ^KE) oc E. It is not generally recognized that the conditions 
kT -C £^1 <C AEf. are quite restrictive and are difficult to satisfy at room temperature. 

The resulting T^^ [Eq. (|^)] under the DK assumptions is proportional to the mobility 
independent of the dominant scattering mechanism [see Eq. (|^)]. In addition, Tf ^ is propor- 
tional to Ef. These trends are not supported by recent experimental measurements ||8|j2T[| 
on 75A n-doped GaAs/Alo.4Gao.6As MQWs at 300K [shown in Fig. |(a,b) (filled circles)]. 
In both cases the experimental trends are weaker than the predicted theoretical ones. Cal- 
culations are shown in Fig. @(a,b) based on our more general theory using OP (solid line) 
and NI (dashed line) as the dominant process determining the mobility. 

Our results agree with experiment if one assumes a shift from OP to NI scattering 
as the mobility drops — this is the origin of the unusual experimental dependence of Ti 
on the mobility. The weaker dependence of Ti on Ei in our theory versus DK theory in 
Fig. 0(b) is due to wavefunction penetration into the barriers and non-perturbative effects. 
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We emphasize as well the key role of temperature studies of the mobility in analyzing the 
temperature dependence of spin coherence. In Fig. §(c) the calculated Ti for one sample 
with a given room-temperature mobility is presented as a function of temperature for NI 
and OP scattering. In particular the OP results appear relatively insensitive to temperature 
from 100-250K — this is due to the rapid decrease in the mobility from OP scattering with 
increasing temperature. This may play a role in the weak temperature dependence seen in 
Ref. 0. 

Figure ^ compares the energy dependence of ^i{E) and ^"^{E) for several additional 
material systems. The cubic dependence of ^1{E) for the three bulk semiconductors is 
confirmed in Fig. |^(a). Fig. ^(b), however, shows that for quantum wells ^KE) is only 
linear (short dashed line for the GaAs MQW) for a small energy range (~ 20 meV) above 
the band edge before it begins to deviate. More energetic states than this certainly contribute 
to the spin coherence times at room temperature. The wider the well the lower the energy 
where ^f{E) deviates from linear behavior, as it approaches a bulk-like E^ behavior. ^1{E) 
for these structures is shown in Fig. 0(c), and is comparable in magnitude to fll{E). As 
the wells become narrower, even the perturbative expressions for Q3 and Qi break down. 
Figure ^(d) shows ^KE) and f2|(-E) for a thin-layer InAs/GaSb superlattice, indicating very 
different behavior from the other structures, poorly reproduced by even the general forms 
of the perturbative expression. 

Table | presents calculations and experimental measurements of Ti for these material 
systems. The order of magnitude discrepancy between DK calculations and measurements 
occurs here as well. Given the uncertainty in experimental mobilities and densities, the 
agreement of our calculations with experiment for both NI and OP scattering is good for all 
systems, and is much better than DK theory. Note that OP and NI scattering calculations 
in the full theory differ from each other by factors of of up to 2 (due to differences in (£")), 
whereas all scattering mechanisms produce the same result in DK theory. As expected, 
for several systems the Ti's are much shorter at higher electron densities, for as the carrier 
distribution is spread further from zone center the effective crystal magnetic fields increase. 
The DK approximation (i) can be evaluated by comparing OPi to OP and NIi to NI, where 
calculations using all terms up to i are designated OP^ and NI^. The difference is up to 
40%. Approximation (ii), however, produces a discrepancy between the DK result and both 
Nil and OPi which can greatly exceed an order of magnitude. 

These calculations consider decoherence arising from the bulk inversion asymmetry (BIA) 
of the constituent materials. We have considered symmetric wells, so another source of 
inversion asymmetry, the structural inversion asymmetry (SIA), does not play a role. In 
other structures, such as single-interface heterostructures, SIA may dominate Interface 
bonding asymmetry (native interface asymmetry, or NIA), which arises in non-common- 
atom structures, could play a role in systems II, IV, or V. The NIA spin splitting for perfect 
interfaces (imperfect interfaces reduce the NIA contribution) has been calculated for System 



II in Ref. |2^. By comparing with Ref. p3| we find the spin splitting of this quantum well 
is dominated by BIA. 

We have presented a quantitatively accurate non-perturbative nanostructure theory for 
electron spin relaxation in bulk and quantum well zincblende semiconductors based on a 
fourteen band model. The calculated electron spin lifetimes in III-V semiconductor bulk 
and quantum well materials are in agreement with experimental measurements, indicating 
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the importance of accurate band structure calculations for zincblende type nanostructures. 

We would like to acknowledge discussions with D. D. Awschalom, T. F. Boggess, 
J. M. Byers, J. M. Kikkawa, and A. Smirl. We would also like to thank D. D. Awschalom and 
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was supported in part by the Office of Naval Research through Grant No. N00014-99- 1-0379. 
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TABLES 



TABLE L Spin coherence times Ti (ps) for several structures, L a 75A GaAs/Alo.4Gao.6As 
MQW i, IL a 70A Ino.53Gao.47As/97A InP MQW [|, IIL an 80A GaSb/80A AlSb MQW 0, 
IV: a 5lA GaAso.19Sbo.8i/8OA AlSb MQW and V: a 2L2A InAs/36.6A GaSb superlattice. 
Calculated times are shown for a given total electron density (n.d. indicates nondegenerate) using 
DK theory (DK), and the nonperturbative theory with optical phonon (OP) and neutral impurity 
(NI) scattering. The subscript i indicates that only terms up to 0,i were used in the calculation. 





System 


Density (cm ^) 


H (cm^/Vs) 


Exp. 


DK 


OPi 


OP 


Nil 


NI 


I 


GaAs/AlGaAs 


2.7 X 10^'^ 


800 


100 


27 


151 


120 


162 


111 


II 


InGaAs/InP 


n.d. 


6700 




1.45 


53 


37 


52 


32 






3.0 X lO^*^ 


6700 


2.6 


0.21 


6.0 


4.9 


6.9 


4.0 


III 


GaSb/AlSb 


n.d. 


2000 




0.59 


1.9 


1.8 


1.5 


1.4 






2.8 X lO^*^ 


2000 


0.52 


0.09 


0.64 


0.55 


0.88 


0.53 


IV 


GaAsSb/AlSb 


n.d. 


2000 




0.09 


0.53 


0.52 


0.44 


0.43 






3.4 X 10^^ 


2000 


0.42 


0.01 


1.9 


1.4 


1.7 


0.87 


V 


InAs/GaSb 


n.d. 


3000 




0.38 


0.77 


0.77 


1.7 


1.6 
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FIGURES 




FIG. 1. T2 in bulk III-V semiconductors as a function of temperature. Solid with squares and 
solid lines respectively represent the results of experiments § and the non-perturbative theory for 
bulk GaAs at the electron density n = 1.0 x 10^^ cm^^. Also shown are results for bulk InAs at 
n = 1.7 X 10^^ cm^^ and bulk GaSb at n = 1.49 x 10^^ cm~^, which are indicated with dashed 
and dot-dashed lines respectively. The difference in slope between GaSb and GaAs occurs because 
GaSb is degenerate for this density. The tabulated mobilities [24| for InAs and GaSb extend only 
to 77K, so at lower temperatures t^{E) was assumed to have the same value as at 77K. 
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FIG. 2. Ti as a function of (a) mobility, (b) confinement energy, and (c) temperature, for 75A 
GaAs/Alo.4Gao.6As MQWs at room temperature. Closed circles represent the results of experi- 
ments 1^. The non-perturbative theory results with OP scattering (solid lines) and NI scattering 
(dashed lines) are shown, as well as the DK theory results (dot-dashed lines). 
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FIG. 3. nj{E) and iljiE) for several structures, (a) ^1{E) for bulk GaAs (solid), InAs 
(dashed), and GaSb (dot-dashed), (b) ^i{E) for GaAs (sohd), InGaAs (long dashed), and GaSb 
(dot-dashed) quantum wells described in Table |. The short-dashed line is the DK approximation 
for the GaAs quantum well, (c) Q'^{E) for the same three structures as (b). (d) J^i(-E) (solid) and 
Q^{E) (dashed) for a thin-layer InAs/GaSb superlattice. 
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